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A shock-like structure appears in a time-averaged pattern produced by the 
Kuramoto-Sivashinsky equation and the noisy Burgers equation with fixed bound- 
ary conditions. We show that the effective viscosity computed from the width 
of the time-averaged shock structure is consistent with that computed from 
the time correlation of the fluctuations. The effective viscosity depends on 
the lengthscale, although our system size is not sufficiently large to satisfy the 
asymptotic dynamic scaling law. We attempt to determine the effective viscosity 
in a finite size system with the projection operator method. 
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1 Introduction 



The Kuramoto-Sivashinsky equation is one of the simplest partial differential 
equations exhibiting spatiotemporal chaos. [jl], ^] In one dimension it has the 
form ^ 

ht — kxx kxxxx ~t~ ~^(hx) j (f) 

where h — h(x,t) is a real function, x G [0,L], and the subscripts stand for 
partial derivatives. An equivalent equation is obtained for u = h x as 

— ^xx Hxxxx ~t~ UUx- (2) 

It is conjectured that the large-scale properties of this equation are described 
by the noisy Burgers equation. ^ || || 

Ut = VkU SC x + \uu sc + £ x (x,t), (3) 

where > represents the effective viscosity on a lengthscale of 1/k (where k 
is the wavenumber), £ x is the spatial derivative of the random noise £(x,t), and 
the parameter A is 1 to satisfy Galilean invariance. The random noise is assumed 
to satisfy (£(x,t)£(x' , t')) — 2Dk5(x — x')8{t — t'). The noisy Burgers equation 
is equivalent to the Karder-Parisi- Zhang equation for a growing interface with 
fluctuations. |(| The dynamic scaling of the noisy Burgers equation (the KPZ 
equation) is expressed as 

(u fc (wK(w)) = k- 3 / 2 g(io/k 3 / 2 ) 

for small wavenumber k and small frequency u. The dynamic scaling of the 
linear Langevin equation 

u t = VkU xx + £ x (x,t) (4) 

for Vk = const and Dk = const is described by (uk(uj)u* k (bj)) — k~ 2 g(u}/k 2 ), 
which differs from the behavior of the noisy Burgers equation. If the effective 
viscosity and the effective noise strength behave as v(k) ~ k^ 1 / 2 , D(k) ~ fc -1 / 2 , 
the power spectrum takes the form k~ z l 2 /{(ui/k 3 ^ 2 ) 2 + I}, which has the same 
form as the dynamic scaling for the noisy Burgers equation. 

The dynamic scaling regime of the noisy Burgers equation has not been 
clearly found in direct numerical simulations of the Kuramoto-Sivashinsky equa- 
tion. Sneppen et al. performed a largescale numerical simulation of the Kuramoto- 
Sivashinsky equation of the form (1).@ They determined the effective diffusion 
constant as v ~ 10.5 and the noise strength as D ~ 6.4 from analysis of the 
interface width ((h(x,t) — (h)) 2 ) in the dynamic scaling regime of the linear 
Langevin equation. They found a crossover toward the dynamic scaling regime 
of the noisy Burger equation. However, they could not demonstarate the dy- 
namic scaling itself. It is believed that a much larger scale numerical simulation 
is necessary to find the dynamic scaling regime of the noisy Burgers equation. 
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Figure 1: Time-averaged shock pattern and a snapshot pattern for the 
Kuramoto-Sivashinsky equation with U = 1.4. 



In this paper, we show that the effective viscosity can be evaluated from the 
relation between the width and amplitude of the time-averaged shock structure 
and the time correlation of the Fourier modes. We do not attempt to find the 
asymptotic dynamic scaling. However, we find some dependence of the effective 
viscosity on the lengthscale. We also performed similar numerical simulations 
for the noisy Burgers equation. We attempt to determine the effective viscosity 
for the Kuramoto-Sivashinsky equation with the projection operator method. 



2 Shock structure in the time-averaged pattern 

We numerically studied the Kuramoto-Sivashinsky equation with fixed bound- 
ary conditions u(0) = —U and u(L) = U in a previous paper, j^] The numerical 
simulations were performed using the Heun method with time step At = 0.002 
and the space step Ax = 0.25. Shock-like structures appear in time-averaged 
patterns for a certain range of boundary values. A nonlinear effect is essential 
for the formation of shock structure. A nonlinear effect appears in this type 
of simulation even for a small system. Contrastly, very large systemsize is nec- 
essary to find a nonliner effect in the analysis of a fluctuating interface width. 
The shock structures are approximated by the functional form 

u(x) = Atanh(K{x - L/2)}. 

Figure 1 displays an example of the shock-like structure for U — 1.4. The 
dashed curve represents a snapshot pattern of u(x, t), the solid curve is a time- 
averaged pattern, and the dotted curve represents 0.58 tanh{0.055(a: — L/2)}. 
It is seen that the amplitude A of the time-averaged shock structure is not the 
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Figure 2: Time-averaged shock patterns and snapshot patterns for the nosiy 
Burgers equation with v = 10, D — 6.5 and (a) U = 3, (b) U = 0.5. 

boundary value U. The effective viscosity v K is evaluated from the relation 
v K = Aj (2k), where A = 1 is assumed, since the space step Ax is sufficiently 
small. (Sneppen et al. estimated A = 2.37 for Ax — 1 and A = 1.025 for 
Ax = 0.5. We believe that the value of A that they found differs from 1, 
because the space step Ax = 1 used in their numerical simulation is too large 
and as a result Galilean invariance is not satisfied.) 

We showed that the effective viscosity depends on the lengthscale of the 
shock width 1/k. The effective viscosity tends to increase as the shock width is 
increased, and it seems to approach a value near 10 in a large lengthscale. The 
value of the effective viscosity is consistent with the value obtained by Sneppen 
et al.§ 

We numerically confirmed that such a time-averaged shock structure appears 
in the noisy Burgers equation (3) with constant bare viscosity Vk — v. The 
boundary conditions were fixed as u(0) = —U and u(L) = U. The numerical 
simulation was performed using the Heun method with time step At = 0.001 
and the space step Ax = 0.25. To maintain u(x)dx = 0, as the case of the 
Kuramoto-Sivashinsky equation, we have further assumed that £ x dx = 0. 
The parameter values are L = 400, v — 10 and D = 6.5. The viscosity and 
the noise strength are chosen to be close to the values found by Sneppen et al. 
Figure 2 (a) displays a snapshot pattern of u(x,t) (dashed curve), the time- 
averaged pattern of u(x,t) = {1/T) J Q u(x,t)dt with T — 50000 (solid curve) 
for U — 3, and a plot of the function tanh{0.15(x — L/2)}. Figure 2 (b) displays 
a snapshot pattern of u(x,t) (dashed curve), the time-averaged pattern (solid 
curve) for U = 0.5, and a plot of the function tanh{0.0227(x — L/2)}. The values 
of it (at) at x — - — - 2.5 x i with integer i are plotted for the snapshot patterns 
to show the three plots clearly. The time-averaged shock structures appear 
clearly also for the noisy Burgers equation, ft is seen that the amplitude of the 
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timc-averagcd shock structure is equal to the boundary value U . The effective 
viscosity can be evaluated from the relation between the shock amplitude and 
the width. The effective viscosities were determined for several values of U, and 
the results are v K = 10 for k = 0.15 (U = 3), v K = 10 for k = 0.05 (U = 1), 
u K = 11 for k = 0.0227 (U = 0.5), and v K = 11.5 for k = 0.013 (U = 0.3). These 
effective viscosities are close to the original viscosity 10. However, they tend to 
increase slightly as k is decreased. 

3 Time correlation of Fourier modes 

It is conjectured that the large scale properties of the Kuramoto-Sivashinsky 
equation are described by the noisy Burgers equation. It is not easy to evaluate 
the effect of the nonlinear term even for the noisy Burgers equation. We study 
the statistical properties of the Kuramoto-Sivashinsky equation by reference to 
the linear Langevin equation (4) with scale-dependent coefficients. The Fourier 
amplitude Uk(t) for the wavenumber k for Eq. (4) obeys 

^ = -v k k 2 u k + f fe (i), (5) 

where £ k (t) is the Fourier transform of £(x,t). The noise £ k (t) is assumed to 
satisfy 

(&(*)) = 0, (£fcWa(0) = 2D k k 2 5(t - t'). 

We further assume that v k and D k may depend on k. 

The stationary probability distribution of u k is Gaussian distribution 

P{u k ) cx exp{-ul/{2a k )}, (6) 

where a k = D k k 2 / '{v{k)k 2 } , and the time correlation function is 

(u k (t)u* k (t')) = <r k exp(-z/ fe fc 2 |i - t'\). 

The decay constant is t = l/{v k k 2 }. We can estimate the effective viscosity 
v k from the decay constant as v(k) — l/(rk 2 ). The power spectrum of u k is 
expressed as (\u k (uj)\ 2 ) ~ D k k 2 /{uj 2 + v k k A }. 

We numerically studied the statistical properties of the temporal fluctuation 
of the Fourier amplitude u k (t) with small wavenumber k for the Kuramoto- 
Sivashinsky equation (2) with periodic boundary conditions. We performed nu- 
merical simulations of the Kuramoto-Sivashinsky equation (2) using the pseudo- 
spectral method with 1024 and 2048 modes. We used system sizes L = 200n 
and 4007T, and the timestep 0.01. Figure 3(a) displays the probability distri- 
bution P{\u k \ 2 ) for k = 0.12, and it is compared with Gaussian distribution 
P(\u k \ 2 ) = l/\/2ira k exp{— \u k \ 2 /(2a k )}, where the numerically obtained value 
of (|«fc| 2 ) is used for a k . The stationary distribution of the temporal fluctuation 
u k can be approximated by the Gaussian distribution. This is consistent with 
the statistical properties of the linear Langevin model. The energy spectrum 
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Figure 3: (a)Probability distribution P(|itfc| 2 ) for k = 0.12 (histogram) and the 
Gaussian approximation Xj^jl-KOk exp{— u 2 k /(2<7k)} (dashed curve) in a semi- 
logarithmic scale, (b) Energy spectrum (|ufe| 2 ) as a function of k. 

(\ u k\ 2 ) ~ Cfc has a peak near k <~ 0.7, which corresponds to the characteristic 
cellular structures in KS turbulence. The energy spectrum is a slightly increas- 
ing function of k in the small- wavenumber region, as shown in Fig. 3(b). This is 
interpreted as meaning that Dk/vk is a slightly increasing function of k in the 
small-wavenumber region. 

The time correlation (uk{t)u* k {t')) was numerically calculated for small k. 
The time correlation function can be approximated by exp(— t/r), except for 
fairly small t. The decay constant Tk was numerically evaluated as the time 
where the normalized correlation (uk(0)u%(t)) / 'ok becomes 1/e. The effective 
viscosity v was evaluated from Vk = l/fak 2 ). Figure 4 displays values of Vk 
as a function of k. We see that the effective viscosity is a decreasing function 
of k. It approaches a value near 10 as k decreases to 0. The effective viscosity 
determined from the width of the shock structure was found to be a decreasing 
function of k, where k is the inverse of the shock width. It approaches a value 
near 10 for small k, and v K is about 5 at k = 0.1. This behavior is consistent with 
the present result, although the lengthscales 1/k and 1/k are slightly different 
quantities. 

We have also calculated the effective viscosity Vk from the time correlation 
for the noisy Burgers equation (3) with v = 10, D = 6.5 and L = 1024 and with 
periodic boundary conditions. Figure 5 displays the numerically evaluated effec- 
tive viscosity Vk- We found that the effective viscosity is almost the same as the 
original viscosity, 10. However, the effective viscosities for small k take slightly 
larger values than the original viscosity. The effective viscosity is expected to in- 
crease as fc~ 1//2 as k — > in the aysmptotic regime of the KPZ dynamic scaling. 
However, this dependence is not clearly seen in this plot. It is expected that 
a much larger scale numerical simulation is necessary to study the asymptotic 
dynamic scaling regime, even for the noisy Burgers equation itself. (Sneppen et 
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Figure 4: Effective viscosity v(k) as a function of k for the Kuramoto- 
Sivashinsky equation. 
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Figure 5: Effective viscosity as a function of k for the noisy Burgers equation. 
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al. determined a critical lengthscale of the crossover from the dynamic scaling 
of the linear Langevin equation to that of the noisy Burgers equation as 2500. 
However, they could not clearly find the aymptotic dynamic scaling regime in 
their numerical simulation of system size 65, 536 for the Kuramoto-Sivashinsky 
equation. It is believed that a very much larger system size may be necessary, 
and we have not yet tried such a larger size simulation.) 

4 Evaluation of the effective viscosity using the 
projection operator method 

We studied statistical properties of the Kuramoto-Sivashinsky equation in the 
previous section. We now attempt to derive the effective viscosity theoreti- 
cally for a finite size system. Recently, Mori and Fujisaka applied a projection 
operator method, which was originally developed for a system near thermal 
equilibrium, to chaotic and turbulent systems. ||] We now apply their method 
to the Kuramoto-Sivashinsky equation. 

The Kuramoto-Sivashinsky equation can be rewritten as 

*£ = (k*-k*)u k +v k , (7) 

where v k — ik/2^2uiu k -i is the nonlinear term. A projection operator P is 
defined as 

Pf({uk}) = (f(u k )u%)/{u k ut)u k , (8) 

where f({u k }) is a function of {ufe}, and the average (• ■ •) is taken with respect 
to the stationary distribution of u k . The operators Q and A are defined as 
Q = 1 - P and 

The operator A is interpreted as the Liouville operator. Applying the projection 
operator method, the Kuramoto-Sivashinsky equation leads to 

dt 



(k z - k 4 )u fc (t) -/ r k (t-s)u k (s)ds + r k (t). (9) 



o 

This equation has the form of a Langevin equation with the memory term 
Jq ^k(t — s)u k (s)ds and the noise term r k {t), where r k (t) and T(s) are expressed 
as 

r k (t) = e tQA Qv k , 

r k (s) = (r k (s)r* k (0))/(u k u* k ). 

The above equations are interpreted as constituting the fluctuation dissipa- 
tion theorem of the second kind. The nonlinearity is included in the equation 
through the terms of T k (s) and r k (t). To evaluate the effective viscosity, we 
have approximated the effective noise term r k {t) as r k (t) = r k (0) +tdr k (0)/dt = 
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Figure 6: Effective viscosity v k as a function of fc for the Kuramoto-Sivashinsky 
equation using the projection operator method. 



Qvk(0) + tQAQv k (0). The effective viscosity of a two-dimensional fluid in ther- 
mal equilibrium was calculated by Iwayama and Okamoto using a similar type 
of expansion. [|lo) 

We have further approximated the integral J Q T k (t— s)uk{s)ds as J ' T k (s)dsu k (t). 
This is a kind of Markov approximation. Then, the effective damping constant 
7^ for u k is approximated as 



7fc 



k 4 



(r k (s)rl(0))ds/(u k u* k ) 



To calculate "f k , we have further assumed that 



r k (s)r* k (0))ds = / {{r k (0)r* k (0))-s(-QAQv k (0)r* k (0))}ds = 



(10) 



(nc(0K(0)> : 



2{-QAQv k (0)r* k (0))' 

where the intergal is approximated at the area of the triangular region satisfy- 
ing (r k {0)r k (0) - s(-QAQv k (0)rl(0)) > 0. The effective damping coefficient 
is expressed as a complicated summation of equal-time correlation functions 
such as (r k (0)r k (0)*) — (Qv k (0)(Qv k (0))*} , although we have used many rough 
approximations. The effective viscosity is evaluated as v k = 7fc//c 2 

We have computed the effective viscosity v k with a numerical simulation, 
where the statistical average (• • •) is replaced by the long time average of the 
chaotic time evolution. Figure 6 displays the numerical results for the system 
size L = 2007T using the pseudo-spectral method with 1024 modes. We obtained 
positive values for the effective viscosity and found that it tends to increase as k 
decreases. However, the effective viscosities found in this way take much larger 
values than the values obtained directly from the time correlation function in 
Fig. 4. This might be due to the crudeness of the approximation. For example, it 
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might be necessary to expand to the next order of t, although more complicated 
summations appear in that case. However, the most essential approximation is 
the Markov approximation, in which it is assumed that the time correlation of 
the noise terms is very short. However, the noise term for the wavenumebr k 
includes the contribution from the slower variables Uk'S with kl < k, and hence 
the approximation of the short time correlation cannot be generally justified. 

5 Summary and discussion 

We have numerically studied the statistical properties of the Kuramoto-Sivashinsky 
equation. We calculated the time-averaged pattern with fixed boundary condi- 
tions and the time correlation of the temporal fluctuation of the Fourier ampli- 
tude with periodic boundary conditions. The effective viscosity v k was evaluated 
from the decay constant of the time correlation function. The numerically ob- 
tained value of Vk is consistent with the effective viscosity v K determined from 
the shock structure, which is due to stepwise fixed boundary conditions. It is 
shown that the effective viscosity depends on the lengthscale. It is interpreted 
as a kind of the fluctuation-dissipation relation for spatio-temporally chaotic 
systems that the response to the external conditions is closely related to the 
temporal fluctuation in the stationary state. Both the response function and 
the temporal fluctuation are described through the effective viscosity. 

The large-scale properties of the Kuramoto-Sivashinsky equation are con- 
jectured to be described by the noisy Burgers equation. We have obtained a 
time-averaged shock structure and the effective viscosity for the noisy Burg- 
ers equation. It is found that there is some correspondence between the two 
equations even in a finite-size system. 

We have attempted to evaluate the effective viscosity for the Kuramoto- 
Sivashinsky equation using the projection operator method. A Langevin-type 
equation was formally obtained using the projection operator method. The 
effective viscosity is formally described in terms of the correlation function of 
noise terms. We have obtained positive values of the effective viscosity. However, 
the values we obtained are larger than the effective viscosity evaluated from the 
time correlation and the shock width. We need to improve the approximation. 

We have not yet found the dynamic scaling regime of the noisy Burgers equa- 
tion in numerical simulations of the Kuramoto-Sivashinsky equation, although 
we have found that the effective viscosity tends to increase as k is decreased. We 
need to study the dynamic scaling regime with numerical simulations of much 
larger systems. 
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